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Abstract 

The issue of how to define and determine an optimal acoustical fit to a set of 
anisotropic elastic constants is addressed. The optimal moduli are defined as those 
which minimize the mean squared difference in the acoustical tensors between the 
given moduli and all possible moduli of a chosen higher material symmetry. The 
solution is shown to be identical to minimizing a Euclidean distance function, or 
equivalently, projecting the tensor of elastic stiffness onto the appropriate symmetry. 
This has implications for how to best select anisotropic constants to acoustically 
model complex materials. 
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1 Introduction 

The theory of wave propagation in anisotropic elastic solids was historically developed 
for homogeneous crystals with a priori known symmetry [1]. The number of elastic 
moduli are then well defined. However, the symmetry is not always apparent in non- 
crystalline and generally complex materials for which measurements can yield as many as 
21 anisotropic elastic coefficients. Whether one has some idea of the underlying material 
symmetry or simply prefers to deal with fewer parameters, the question arises of how to 
best fit the given set of elastic moduli to, for instance, a transversely isotropic material 
model. This issue occurs in acoustical measurements of composites composites [2] , and in 
geophysical applications [3] where laboratory measurements might yield 21 moduli, but 
seismic modeling requires a higher symmetry, such as transverse isotropy. The purpose 
of this paper is to provide a simple but unambiguous means to find the reduced set of 
anisotropic elastic constants that are in a certain sense the best acoustical fit to the 
given moduli. The optimal material minimizes the mean square difference in the slowness 
surfaces of the given moduli and of all possible sets of elastic constants of the chosen 
symmetry. 

The prevailing approach to finding a reduced set of moduli does not invoke acoustical 
properties, but views the moduli as elements in a vector space which are projected onto 
the higher elastic material symmetry. This is achieved by defining a Euclidean norm 
for the moduli C according to || C|| 2 = CijuCijki where 1 C^-jy are the elements of the 
stiffness tensor. This provides a natural definition for distance, from which one can 
find the elastic tensor of a given symmetry nearest to an anisotropic elastic tensor, or 
equivalently, define a projection appropriate to the higher symmetry. Gazis et al. [4] 
outline the procedure in terms of fourth order tensors, while more recently Browaeys and 
Chevrot [3] provide projection matrices for C expressed as a 21-dimensional vector. Helbig 
[5] provides a useful overview of the problem from a geophysical perspective, Cavallini [6] 
examines isotropic projection specifically, and Gangi [7] gives formulas for several other 
symmetries. One drawback of the approach is that while it minimizes the Euclidean 
distance between the original and projected stiffness tensors, it does not provide the 
analogous closest compliance (the inverse of stiffness). In this sense the stiffness projection 
method, while simple and attractive, is unsatisfactory because it is not invariant under 
inversion. Alternative procedures based on non-Euclidean norms such as the Riemannian 
[8] or log-Euclidean [9] metrics do not have this deficiency, and, in principle, provide a 

1 Lower case Latin suffices take on the values 1, 2, and 3, and the summation convention on repeated 
indices is assumed. 
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unique "projection" regardless of whether one uses the stiffness or the compliance tensor. 

An apparently quite different approach is to try to find higher symmetry moduli which 
in some way better approximate the acoustic properties of the given moduli. Thus, Fe- 
dorov [10] considered the question of what elastically isotropic material is the best acoustic 
fit to a given set of anisotropic moduli. He defined best fit to mean the effective bulk 
and shear moduli {k, fi} which minimize the mean square difference between the slowness 
surfaces of the original anisotropic material and the isotropic material characterized by 
{k, /i} (density is unaffected). Fedorov obtained explicit expressions for the moduli, eq. 
(26.19) of [10], or in the present notation 



Fedorov's procedure for finding a suitable set of higher symmetry moduli is physically 
appealing, especially as it seeks to approximate acoustical properties. Also, as Fedorov 
and others [10, 4, 3] have shown, k and fi are precisely the isotropic moduli found by 
the stiffness projection method. However, Fedorov only considered effective isotropic 
moduli and it does not appear that anyone has attempted to generalize his method to 
symmetries other than isotropic. The purpose of this paper is to solve what may be 
termed the generalized Fedorov problem. The solution is found for the acoustically best 
fitting moduli of arbitrary symmetry, which is of higher symmetry than the given moduli 
but lower than isotropy. The central result is that the solution of the generalized Fedorov 
problem possesses the same crucial property as the solution obtained by Fedorov [10], 
that is, the best fit moduli are identical to those obtained by the stiffness projection 
method. This result provides a strong physical and acoustical basis for using the Euclidean 
projection scheme that has been absent until now. 

We begin in Section 2 with the definition of elastic tensors and associated notation. 
The generalized Fedorov problem is introduced and solved in Section 3. Examples are 
given in Section 4. 

2 Preliminaries 

2.1 The elasticity tensor and related notation 

The solution of the generalized Fedorov problem is most easily accomplished using ten- 
sors, which are reviewed here along with some relevant notation. Boldface lower case Latin 
quantities indicate 3-dimensional vectors, such as the orthonormal basis {e^, i = 1,2,3}, 
"ghostface" symbols such as C indicate fourth order elasticity tensors, and boldface 
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capitals, e.g. A, are second order symmetric tensors, with some exceptions. Com- 
ponents defined relative to the basis vectors allow us to represent arbitrary tensors in 
terms of the fundamental tensors formed from the basis vectors, thus, A = ® e^, 

C = Cijki e,i (g) e.j ®e.]~® e^, where ® is the tensor product. The identity tensors of second 
and fourth order are I and I, respectively, and we will also use the fourth order tensor 
J. These have components 

hi = hjki = -^{&ik&ji + 8u5jk), JijM = -^8ij5ki ■ (2) 

Isotropic elasticity tensors can be expressed as a linear combination of I and J. Using 
Lame moduli /i and A = k — -//, for instance, 

C = 2/iI + 3AJ. (3) 

Products of tensors are defined by summation over pairs of indices: ( CA)jj = C^iA^i 
and {AM) ijM = A ijpq B pqkl . Thus, II — JI = I, II = I, J J = J, and I J = J I = J. 
The inner product of a pair of tensors of the same order is defined 

(u, v) = tr(uv), (4) 

where tr A = An and tr A = A^ , and the Euclidean norm of a tensor is 

\\u\\ = (u, u) 1 / 2 . (5) 

This will be used to compare differences between tensors. 

An anisotropic elastic stiffness tensor C relates stress T and strain E according to the 
generalized form of Hooke's law: 

T = CE . (6) 

Stress and strain are symmetric second order tensors, implying that their components are 
also symmetric, = Tji, = Eji, which in turn implies the first two of the following 
identities 

Cijki — Cijik, Cijki = Cjiki, Cijki = Ckuj- (7) 

The third property is a consequence of the existence of a positive strain energy of the form 
W = |(E, CE), which also constrains the moduli to be positive definite. The squared 
norm of the elasticity tensor is 

|| C|| 2 = ( C, C) = CijkiCijki ■ (8) 
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The 21-dimensional space of elasticity tensors can be decomposed as a 15-dimensional 
space of totally symmetric fourth order tensors plus a 6-dimensional space of asymmetric 
fourth order tensors [11, 12], by 

C = C (s) + C (a) , (9) 

where 

(M. = 

3 



Qjkl — a (Pijkl + Cukj + Cikji). (10) 



The elements of the totally symmetric part satisfy the relations C\ s ^ kl = C^ 7 in addition 
to (7). Thus, C^li is unchanged under any rearrangement of its indices, and has 
at most 15 independent elements. This can be seen by the explicit representation of 
the asymmetric part in terms of the remaining six independent elements, which are the 
components of the symmetric second order tensor 

Dij = Cijkk ~ Cikjk- (11) 

Thus, 

Cijli = g [ZDijlki + 2Dkihj — Diklji — Duljk — Djkhi — Djilik + D mm (Iijki — 37^)] . (12) 

An asymmetric tensor satisfies, by definition [11], 

C£i + C& + Cgi = <). (13) 

Note that the symmetric part of an asymmetric tensor is zero, and vice versa, and that the 
decomposition (9) is orthogonal in the sense of the Euclidean norm, || C|| 2 = || C (s) || 2 + 
|| C (a) || 2 . The partition of C as a sum of totally symmetric and traceless asymmetric 
tensors is the first step in Backus' harmonic decomposition of elasticity tensors [11], a 
partition that has proved useful for representations of elastic tensors [13] and has also 
been used to prove that there are exactly eight distinct elastic symmetries [14]. Finally, 
for future reference, note that the totally symmetric part of the fourth order identity is 

lifki = l( S K S ki + SikS j i + S a S kj ) & I« = ^I+J. (14) 
2.2 The acoustical tensor and the tensor C* 

The acoustical tensor, also known as Christoffel's matrix [1], arises in the study of plane 
crested waves with displacement of the form u(x, t) — b h(n ■ x — vt), where b is a fixed 
unit vector describing the polarization, n is the phase direction, also a unit vector, v is 
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the phase velocity, and h is an arbitrary but sufficiently smooth function. Substituting 
this wave form into the equations of motion 

d 2 u k .. 

c - ki d^r pUi ' (15) 

where p is the density, implies that it is a solution if and only if b and pv 2 are eigenvector 
and eigenvalue of the second order tensor Q, 

Q ik (n) = Cijiarij-ni . (16) 

This definition of the acoustical tensor is not the product a fourth order tensor with a 
second order tensor. In order to express it in this form, which simplifies the analysis later, 
introduce the related fourth order tensor C* defined by 

Qjki = 2 {C ik ji + Cujk) • (17) 

Thus, 

Qij(n) = C* jkl n k ni <=> Q(n) = C* n ® n . (18) 

The operation defined by * is of fundamental importance in solving the generalized 
Fedorov problem, and therefore some properties are noted. First, * is a linear operator 
that commutes with taking the symmetric and asymmetric parts of a tensor: C^-** = 
C *W = C («) and c <«0* = C *(o) = _i C W. Accordingly, C* is partitioned as 

C* = C {s) - ^C (a) , (19) 

and repeating the * operation n times yields C"* = + (—1/2)™ C^ a \ Taking n = 2 
implies the identity 

C = 2C**-C*, (20) 

from which C can be found from C*. Hence the mapping C <-> C* is one-to-one and 
invertible, that is bijective. This property is important in the inverse problem of determin- 
ing elastic moduli from acoustic data [15]. Acoustic wave speeds and associated quantities 
are related primarily to C* through the acoustical tensor also known as the Christoffel 
matrix [1], and this can be determined uniquely from C using (20). Decompositions of C 
and C* into totally symmetric and asymmetric parts are unique, and knowledge of one 
decomposition implies the other. In particular, C = C* if and only if the asymmetric 
parts of both are zero. This occurs if the moduli satisfy CV,- fcfc = c7j fcjfc , which together 
with the symmetries (7) are equivalent to the Cauchy relations, see Sec. 4.5 of Musgrave 
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[1]. We note that the operation * is self adjoint in the sense that the following is true for 
any pair of elasticity tensors, 

(A, B*) = (A*, B) . (21) 

Elastic moduli are usually defined by the Voigt notation: Cijki = cjj, where /, J = 
1, 2, 3, ... 6 and 1—1,2, 3, 4, 5, 6 correspond to ij = 11, 22, 33, 23, 13, 12, respectively, i.e. 



C 



c 



feu 


Cl2 


Cl3 


C14 


Cl5 


ci 6 \ 






C23 


C24 


C25 


c 26 






C33 


C34 


C35 


C36 








C44 


C45 


C46 


S 


Y 


M 




C55 


C56 


{ 










c 66/ 



(22) 



The components of C* are then 



C* 



Cll 


c 66 


c 55 


c 56 


Cl 5 


Cl6 




C22 


C44 


C24 


C46 


C26 






C33 


C34 


C35 


C45 



s 



Y 



M 



^(044 + 023) |(C45 + C36) 
|(c 5 5 + C13) 



|(C46 + C25) 
|(C 5 6 + C14) 
|(c 66 + C12)/ 



(23) 



3 Fedorov's problem for particular symmetries 
3.1 Definition of the problem 

We assume as known the elasticity tensor C of arbitrary material symmetry with as many 
as 21 independent components. A particular material symmetry is chosen with prescribed 
symmetry axes or planes. For instance, transverse isotropy with symmetry axis in the 
direction a, or cubic symmetry with orthogonal cube axes a, b, c. Fedorov's problem for 
particular symmetries is to determine the elastic stiffness Cs ym of the chosen material 
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symmetry which is the best fit in the sense that it minimizes the orientation averaged 
squared difference of the acoustical tensors. We introduce the acoustical distance function 

2-7T 7T 

f(C,C Sym ) = ^- J dcf> J wx9d9\\Q(C,n) - Q(C Sym ,n)\\ 2 , (24) 
o o 

with n = sin #(cos 0ei + sin <pe2) + cos 9e^. The same function was considered by Fedorov 
[10] with Csym restricted to isotropic elasticity Thus, substituting Cs ym = «i I+«2 J into 
(24) one obtains, after simplification, a positive definite quadratic in the two unknowns 
a\ and a?,- A simple minimization yields a± = 2/i and «2 = 3k — 2/i where k and fj, are 
defined in eq. (1). The symmetry of Cs y m can be considered as arbitrary for the general 
problem addressed here, although Fedorov's isotropic result is recovered as a special case 
of the general solution discussed in Section 4. 

The integral over directions in (24) can be removed using 

2-7T 7T 

— / d</> / sin 6d9 n®n®n®n = - (25) 
o o 

where 1^ is defined in eq. (14). The identity follows by noting that the integral must 
be a totally symmetric isotropic fourth order tensor of the form al^ s \ Taking the trace 
of both sides and using tr 1^ = I^Jj = 5, tr(n (gin ® n ® n) — 1, gives a = 1/5. Thus, 
since 

Q(C,n) = C*n®n, Q( C Sym , n) = C* Sym n ® n, (26) 
the distance function reduces to 

/ = ~ ((C* - C^J, lW(C* - C^J) . (27) 

Define the modified inner product for elasticity tensors: 

( A, B)„ = tr( I W AB) = ( A, B) , (28) 

and norm 

||A|| a = (A, A)y 2 . (29) 
Then Fedorov's problem for particular symmetries amounts to: 

Fedorov <^ minimize || C* — C* Sym \\ a ■ (30) 



The reason the problem is expressed in this form is to make the connection with the 
notion of projection onto the chosen elastic symmetry, or equivalently of finding the elastic 
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tensor of the chosen symmetry nearest to the given elasticity C. As mentioned in the 
Introduction, this question has been addressed by several authors and has an explicit 
solution. The issue is to find the elastic tensor Cs ym which minimizes the Euclidean 
distance function d, where 

d{ C, C S ym) = || C - C S ym\\ ■ (31) 

Comparing eqs. (30) and (31), the two problems appear tantalizingly similar. It should 
be realized that Fedorov's problem involves C* not C, and that the norms are distinct. 
However, it will be proved that the problems share the same solution: 

{ C S ym : min || C* - Cj, ym || a } = { C Sym ■ min || C - C Sy m\\} , (32) 

which gives the central result of this paper, i.e. 

Fedorov <^ Euclidean projection. (33) 

This equivalence enables us to provide an explicit solution to the generalized Fedorov 
problem, e.g. using the methods of Gazis et al. [4], Browaeys and Chevrot [3] or others. 
The remainder of this Section develops a proof by construction of the solution of the 
Fedorov problem, which is shown to be identical to the Euclidean projection. Some 
further concepts and notation are required and introduced next. 

3.2 Basis tensors 

2 The solution uses a fundamental decomposition of the chosen material symmetry us- 
ing basis tensors. These form a vector space for the symmetry in the sense that any 
elasticity tensor of that symmetry may be expressed uniquely in terms of N linearly in- 
dependent tensors Vi, V2, • • ■ Vat, where 2 < iV < 13 is the dimension of the vector 
space for the material symmetry. For instance, isotropic elasticity tensors are of the form 
Ciso = «i 1 + a 2 J, «i, a<i > 0. The procedure is analogous for other material symmetries, 
cubic, transversely isotropic, etc., and is described in detail by Walpole [16] who pro- 
vides expressions for the base tensors of the various symmetries. Thus, N = 2, 3, 5, 9 for 
isotropy, cubic symmetry, transverse isotropy and orthorhombic symmetry, respectively. 
N = 13 corresponds to monoclinic, which is the lowest symmetry apart from triclinic 
(technically N — 21) which is no symmetry. The precise form of the basis tensors is 
irrelevant here (examples are given in Section 4, full details are in [16], and Kunin [17] 
2 
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develops a similar tensorial decomposition), all that is required is that they be linearly 
independent, and consequently any tensor with the desired symmetry can be written 

N 

c Sym = J2& Y *- ( 34 ) 

i=X 

The coefficients follow by taking inner products with the basis tensors. Let A be the the 
N x N symmetric matrix with elements 

A i3 - = (V i5 V,). (35) 

A is invertible by virtue of the linear independence of the basis tensors, and therefore 

N 

Conversely, eq. (34) describes all possible tensors with the given symmetry, and in par- 
ticular, Csym = if and only if fa — 0, i — 1, 2, . . . N. 

3.3 The elastic projection 

It helps to first derive the solution that minimizes the Euclidean distance function of 
(31). By expressing the unknown projected solution in the form (34), it follows that the 
minimum of || C — Cs^H 2 is determined by setting to zero the derivatives with respect 
to Pi, which gives the system of simultaneous equations 

]T(V t , = C), i = l,2,. ..N. (37) 

The inner products are the elements of the invertible matrix A, and so 

N 

C S ym=Y2(C, VjA^Yj. (38) 

Furthermore, the distance function at the optimal Cs ym satisfies 

||C-C 5 , m || 2 =||C|| 2 -||C 5 , m || 2 , (39) 

as expected for a projection using the Euclidean norm. This is essentially the method 
used by Arts et al. [18]. Helbig [5] describes the procedure as the optimum approximation 
of an arbitrary elasticity tensor by projection of a 21- dimensional vector onto a subspace 
of fewer dimensions, and Browaeys and Chevrot [3] list the explicit forms of the 21D 
projection operators for various symmetries. 
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Define the projection operator Vsym which maps C onto the chosen symmetry, 

TV 

V Sym C=J2(C, VJAjVj. (40) 

Equations (38) through (40) imply that Vs ym C is the Euclidean projection, also equal to 
the closest elasticity tensor of the chosen symmetry to C. We note the following important 
property: 

V Sym C* = {VsymCy. (41) 

In other words, the operation * commutes with the projection operator. This is not 
surprising if one considers that the * operation is a linear mapping on the symmetric and 
asymmetric parts of of C, and therefore * the maintains the material symmetry of C. 
However, a more detailed proof of the identity (41) is provided in the Appendix. 

3.4 Solution of the generalized Fedorov problem 

We now calculate the optimal Cs ym for the generalized Fedorov problem, and show that 
it is equivalent to the moduli from the Euclidean norm. The starting point this time is 
to express the unknown C* Sym (rather than Cs ym ) in terms of the basis tensors, 

TV 

C^, = E«.V, (42) 

i=l 

This is justified by the fact that all tensors of the given symmetry are linear combinations 
of the basis tensors. Furthermore, the coefficients on are related to those in eq. (34) by 
the N x N matrix P introduced in the Appendix. Let a and (3 denote the N x 1 arrays 
with elements ctij and a iy then 

(3 = P*^ <=> a = 2P t (3 - (3. (43) 

Equation (30) implies that the coefficients coefficients ai, a.2, ■ ■ - oln satisfy the system of 
iV equations 

N 

^(V. t , Wj)aaj = {V u C*) a , i = l,2,...N, (44) 
or, in matrix format, 

SAa = S 7 . (45) 
Here 7 with elements 7« and the N x N matrix S are defined by 

7* = (V*, C*>, (46a) 

N 1 

% V, = - ( I {s) V, + V, I (s) ) = U< . (46b) 
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The tensors U$, i — 1, 2, . . . N form a linearly independent set of basis tensors for the 
given symmetry. This may be seen by assuming the contrary, i.e. that there is a set of 
non-zero coefficients a, such that 

TV 

^a 4 U,=0. (47) 

i=l 

Let A be the non-zero tensor 

N 



5>V,-, (48) 



i=l 



then eq. (47) requires that 



-A+~(AI)<g)I+~I®(Al r ) = 0. (49) 
Multiplication by J implies 

1 3 3 

-AI+-(J, A)/ = 0, (50) 
o 2 

and taking the inner product with i" gives 

11 (J, A) = 0. (51) 

Therefore, 

A = 0, (52) 

and the U, tensors form a linearly independent basis for the symmetry. In particular, 
the tensors Vj can be expressed in terms of this alternate basis, and so the matrix S is 
invertible. Hence 

a = A" 1 7. (53) 

We are now in a position to determine the optimal C* Sym and hence Cs ym - Equations 
(42) and (53), along with (40), imply 

N 

C *Sym = E< C *' ViJAT. 1 ^ 

= VsyrrX*- (54) 

Using the fundamental property of the projection operator (41), the optimal elasticity as 
determined by (54) is seen to be exactly the same as the Euclidean projection, i.e. of eq. 
(38). This completes the proof of the main result, the equivalence (33). 



4 Examples 

The general procedure for projection is illustrated for several symmetries, and an example 
application is discussed in this Section. 
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4.1 Basis functions for isotropic, cubic and hexagonal materials 

The tensor decomposition procedure is described for the three highest symmetries: isotropic, 
cubic and transversely isotropic. The fundamental matrices A and P are given explicitly. 

4.1.1 Isotropic approximation (N = 2) 

Let the basis tensors be Vi = J, V2 = K = I — J. Then A is a 2 x 2 diagonal matrix, 
A = diag(l, 5) and the optimal moduli are 

Ciso = 3k J + 2/i K, (55) 

where 

«=i(J, C), /i=l(K, C). (56) 

This is precisely Fedorov's original result [10], eq. (1). The * operation is defined by the 
matrix P of eq. (Al) 1; which is 



1 1 

3 3 

5 1 

3 6 



(57) 



4.1.2 Cubic materials (N = 3) 

Let a, b, c be the cube axes, and select as basis tensors [16] Vi = J, V2 = L = I — H, 
and V 3 = M = H - J, where 

H = a®a<g)a®a + b<g)b®b<g)b + c®c®c®c. (58) 

Then A = diag(l, 3, 2) and the optimal moduli of cubic symmetry are 

C Cub = 3k' J + 2// L + 27] M, (59) 

where 

K?=UlC), n' — ^(L, C) V = ](M,C). (60) 
3 6 4 

It is assumed here that the axes of the cubic material are known, otherwise a numerical 
search must be performed to find the axes which give the closest, i.e. largest, projection. 
This additional step is discussed in the numerical example below. Also, 



/ I I I \ 
'3 3 3 » 



1 _ 1 

2 2 



.2 _i 2 

\ 3 33 



(61) 
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4.1.3 Transverse isotropy /hexagonal symmetry (N = 5) 

Let a be the direction of the symmetry axis, and define [16] 

Y 1 = A®A, Y 2 = ±B®B, V 3 = \{A®B+B®A), V 4 = 2V*, V 5 = 2V;-V 2 , 

(62) 

where A = a (g) a and B = I — A. Then A = diag(l, 1, 1, 2, 2) and the optimal TI 
moduli are given by 



3 ^5 



(63) 



i=4 



Also, 



/ 1 \ 



| \ 



\ 



i I 



(64) 



\ 1 / 
4.2 Application to acoustically measured data 

Let us assume that the material has cubic symmetry but the cube axes orientations are 
unknown. The effective cubic moduli are, in coordinates coincident with the cube axes, 
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Cub 



Ml 


C 12 


C 12 











L 12 


c c u 


C 12 











r c 

c 12 




L ll 




















C 66 




















C 66 



















C 66 



(65) 



where the three constants can be expressed in terms of the bulk modulus and the two 
shear moduli, 

4 2 
c c u = «' + 3^ c u = «' - 4e = • ( 66 ) 
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The effective bulk modulus of (60) is an isotropic invariant that is independent of the 
cube axes orientation, 

K ' = \ C Hjj = g ( C H + C 22 + C 3 3 + 2c 12 + 2c 23 + 2c 31 ) . (67) 

Similarly, the combination (3/i' + 2r])/5 = /i, the isotropic shear modulus, implying 
1 2 

6/i' + 4?7 = C m - -C iij:j = - (en + c 22 + c 33 - C12 - c 23 - c 3 i) + 2 (c 44 + c 55 + c 66 ) . (68) 

Only one of the three cubic parameters depends upon the orientation of the axes, and it 
follows by considering the squared length of the projected elastic tensor, 

II CcubW 2 = 9k' 2 + 12/i' 2 + 8r/ 2 = 9k' 2 + -(6// + 4r/) 2 + -(2// - 2r/) 2 . (69) 

5 5 

The first two terms are independent of the cube axes orientation, and therefore the closest 
cubic projection maximizes the final term. Note that 

2/i' -2 V = \c m + -C lin - \{ C, H) = \k' + 2/i' + \r) - \{ C, M> . (70) 
o 6 6 2 d 6 

The moduli k', fi f and r] must all be positive in order for the material to have positive 
definite strain energy. Therefore, (2/i' — 2rf) 2 is maximum when ( C, H) is minimum. The 
latter is also a positive quantity, which may be expressed 

( C, H) = Ci jk iHi jk i = c' n + c' 22 + c 33 , (71) 

where c'jj are the elements of C of (74) in the coordinate system coincident with the cube 
axes. 

We can also write 

(C, H) = tfCa + tfCb + tfCc, (72) 

where 

a! = [a\, a 2 ,, a 3 , 2a 2 a 3 , 2a 3 ai, 2aia 2 ), etc. (73) 

The minimum can be found by numerically searching over all possible orientations using 
Euler angles to transform from the fixed vector basis to the cube axes. The rotated (cube) 
axes are simply the columns of the 3x3 transformation matrix R G 5*0(3). 

As a specific application we reconsider the elastic moduli obtained from ultrasonic 
measurements by Francois et al. [2]. In the notation of eq. (22), the raw data for the 
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stiffness tensor is 3 
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(74) 



We find that 



C 



Cub 



/213.4 148.5 148.5 

148.5 213.4 148.5 

148.5 148.5 213.4 







\ 






139.8 
139.8 



(GPa), 



(75) 



\ 139.8/ 

in agreement with Francois et al. [2]. The rotation from the coordinate system {ei, e2, 63} 
of (74) can be expressed as a rotation about a single axis n through angle 6 using Euler's 
theorem [19]. The angle and axis follow from 2 cos 6* = tr R — 1 and 2 sin# n = skew R, 
where skew Y = —eijkYij and is the third order alternating tensor. Thus, we find 
9 = 96° and n = (0.18, 0.06, 0.98), and the full set of moduli in the rotated frame are 

/228 141 148 -5 -1 -2\ 



a 



141 209 156 -6 23 -1 

148 156 203 -7 -2 4 

-5 -6 -7 144 12 -3 

-1 23 -2 12 139 11 

\-2 -1 4 -3 11 136/ 



(GPa). 



(76) 



3 The element C53 is given as 49 in [2], which appears to be a typographical error. 
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This set of cube axes is unique within rotation under the group of transformations con- 
gruent with cubic symmetry. In this preferred coordinate system, the projection onto the 
cubic moduli in this frame is simply 

4 = K c 'u + 4 + 4) =213.4 (GPa), 

4 = 1(4 + 4 + 4) =148.5, (77) 
C 66 = |( C 44 + c 55 + c m) = 139.8. 

5 Conclusion 

The main result of this paper is the proof that the Euclidean projection of anisotropic 
elastic constants onto a higher material is identical to minimizing the mean square dif- 
ference of the slowness surfaces. This provides a well grounded acoustical basis for using 
the Euclidean projection as a natural way to simplify ultrasonic or acoustic data. The 
equivalence generalizes the original result of Fedorov for the best isotropic acoustical fit 
to a given set of anisotropic moduli. 
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Appendix 
A. The matrix P and proof of eq. (41) 

The tensor V* can be expressed using eq. (19) in terms of the totally symmetric and 
asymmetric parts of Vj. Both and possess the same material symmetry as Vt, 
and hence the symmetry is inherited by V*. Since the Vj themselves form a basis for the 
material symmetry, it follows that V* can be written as a linear combination of them. 
Let P be the N x N matrix which defines the * operation in terms of the basis tensors, 
that is, 

N 

v* ^/V4 (a.i) 

It follows from eq. (20) that 

N 

V l = -V* + 2^P JJ V*, (A.2) 
j'=i 
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and hence the inverse of P is 

P- 1 = 2P-I W , (A.3) 

where I(jv) is the N x N identity. 
Consider the identity 

(V*, V i > = <V i) V*), (A.4) 

which follows from (21). Using (Al)i to eliminate V* and V* from the left and right 
members, eq. (A.4) implies 

PA = AP*. (A.5) 

It is worth noting the very specific nature of the matrix P that is required to satisfy eq. 
(A.5). This matrix essentially defines the * operator, from which the totally symmetric 
and asymmetric parts of a tensor can be found in terms of the basis tensors. Some 
examples of P are given in Section 4. 

We now turn to the proof of eq. (41). Starting with the definition of Vs ym in (40), we 
have 

N 

»,3=1 
N 

N 

= ]T(C, YjXijY*, (A.6) 

where (21) has been used in the second line, and the matrix X is defined as 

X = P'A~ 1 P- 1 . (A.7) 

Comparison with (Al) 2 gives 

X = A" 1 , (A.8) 

and substituting from (A.8) into (A.6) and comparing it with (38) implies the fundamental 
property (41). 
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